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FOREWORD 
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personnel in the Fluid Mechanic Applications Group of the Lockheed- 
Huntsville Research & Engineering Center for NASA-Marshall Space 
Flight Center. 
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(M553) and R. M. Poorman (M551) of NASA-MSFC S&E-PE-MW and 
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SECTION I. INTRODUCTION 


Lockheed -Huntsville began work in May 1972 to assist the 
National Aeronautics and Space Administration, Marshall Space 
Flight Center, in meeting the objectives of the M551 Metals Melting 
Expo'-iment and the M553 Sphere Forming Experiment to be conducted 
in the M512 Facility aboard the Skylab Laboratory. Specifically, 
ground-based studies have been conducted with the objectives of 
clarifying the effects of gravity and optimizing the experiments within 
the constraints of existing hardware. All process phenomena have 
been considered with particular attention given to adhesion-cohesion 
studies with emphasis on the fluid dynamics of the mo’ten metal. 
Detailed thermal histories have been generated for the three- 
dimensional specimen geometry including radiation, conduction and 
vaporiration losses with allowances for variable properties. In 
addition, quantitative tecliniques have been used to establish the 
gravitational level (magnitude and direction) for KC-135 aircraft 
flights and the actual Skylab mission for both experiments. 

The overall study is composed of three distinct portions. 

Phase A — the preparation of a ground base study plan. Phase B — 
the performance of the study plan, and Phase C — the analysis of 
experimental flight data. This final report documents the results 
of the overall effort. 

The M551 Metals Melting Experiment consisted of a rotating 
metal disk mounted perpendicular to an electron beam (eb) heat source 
with the eb impingement point located 6 cm from the center of rotation. 
The velocity of the beam relative to the impingement point was 1.61 
cm/sec with a corresponding radial acceleration of 5.09(10)“^g. The 
eb was focused to approximately 0.15 cm diameter. The three mate- 
rials listed below were used in the experiment. 

• 2219 Aluminum 

• 321 Stainless Steel 

• Tantalum. 

The thickness of the material varies with angular position around 
the disk. 

The M553 experiment consisted of two 15 -specimen indexing 
pinwheels mounted perpendicular to the electron beam. With each 
discrete rotation of the pinwheel a new cylindrical specimen was 
positioned in the path of the electron beam. As the specimen melted, 
the hardware was designed to allov/ for deployment of the resulting 
liquid spheres with the ensuing free-float time giving rise to a con- 
tainerless solidification in the vacuum chamber. 
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Before presenting specific calculations and experimental results 
relating to the M5 51 and M553 experiments, a brief discussion of 
some fundamental concepts will be presented (Section II). Specific 
attention will be given the M551 experiment in Section III, while 
Section IV is devoted to the M553 experiment. 
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SECTION II. FUNDAMENTAL CONCEPTS 


This section reviews some of the fundamental concepts pertinent 
to the M551 and M553 experiments. Specific calculations relative to 
either the M551 or M553 experiment are presented in Sections III or 
rv, respectively. 


A. MECHANISMS OF POWER LOSS 

The distribution of eb power can be described by the following 
equation: 


Loss 


L + 9c «r + % 


( 1 ) 


where 


‘i-/ 

‘ir 

•ii 


power loss 

power loss 

power loss 

power loss 
evaporated 


due to vaporization 
due to heat conduction 
due to radiation 

due to ionization and excitation of 
atoms. 


The vaporization loss can be calculated from the Langmuir 
theory. Here, the mass flux due to vaporization can be calculated 
from kinetic theory and is equal to 


rh = P 

V 



(2) 


where 


P = vapor pressure 
M = molecular weight 
R = universal gas constant 
T = temperature. 
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If more than one component is present, Raoult's law may be 
used to calculate the equilibrium vapor pressure that is exerted by 
a component in a solution. This law states that 


^a 


N P 
a V 

a 


n 


, a D c 




( 3 ) 


where 


n , n, , n 
a’ b’ c 

N 

a 

1 he 


= partial pressure components "a" in a solution 
wit h, b, c, . . 

- vapor pressure of "a" in pvire state 

= moles for components a, b, c 
= mole fraction of "a." 

total vaporization power loss can now be calculated from 



( 4 ) 


where A is the energy required to vaporize a unit mass at 0°K (i.e., 
the sum of the heat of melting and the heat of vaporization). 


A second mechanism of heat loss is due to thermal conduction. 
This can be obtained (via a thermal analyzer) from 


q^ =y* q • n dA (5) 

A 

where q is the heat flow vector and n is the surface normal. The 
term q is calculated by 


q = -kVT 


with T satisfying 


V 


i AX 

a at 


» 

Q 

k 


( 6 ) 


( 7 ) 
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where 

a = thermal diffusivity 
k = thermal conductivity 
Q = heat generation (source). 

The heat loss to radiation can be calculated via the Stefan- 
Boltzmann law 


-f 


4 

e a T y dA 


(8) 


where 

e = emissivity 

y = fraction of radiation differential area dA which 
escapes from eb hole 

a = Stefan-Boltzmann constant. 

The average power lost to ionization and excitation can be 
calculated from the universal electron range (R) — incident electron 
energy relationship. This can be written as 

( 9 ) 
( 10 ) 


dE 

dR 


= known constant 


The density of the gas in the hole is 


m 

P = IT- 

^o^ 


whe. •- 


A^ = area of cavity 
V = average molecular velocity. 
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The energy loss due to ionization and excitation is then 

= ( 11 ) 

whe re 

I = current 
d = cavity depth. 

These mechanisms are shown schematically in Figure 1. 

D. ACTIVE FORCES ATTRIBUTABLE TO ELECTRON BEAM 
Electron Force 

The transfer of electron momentum during impingement on a 
surface can be eqxiated to the net force on the surface and is given by 


F /a ^ 
e 


J 

m V — 


e e e 


where 


F /a 
e' 


m 

e 


e 

e 


J 

The velocity 


= force per unit area 
= mass of electron 
= velocity of electron 
= charge of electron 
= current density 
can be calculated from 


V 

e 



where 

E = acceleration potential. 


( 12 ) 


(13) 
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Thus llie total force on area A is 


TmE 

e 


1 7 

with e/m = 5.27(10) slat coulomb/gm. 


(14) 


Vaporization Force 

The reaction force of evaporating mass can be calculated from 
conservation of momentum: 


F /a = m V 

V 


wliere 


(15) 


V = molecular velocity. 

The most probable velocity calculated from a Maxwell-Boltzmann 
distribution is 



while the mean speed is given by 



where 


k = Boltzmann's constant 
m = molecular mass 


The total force can be found by integrating over the surface 


=y'rh (T) V(T) dA. 
A 


(18) 
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T 


a 


If the most probable velocity is used this equation reduces to 

*'v = rk/ Pv'T* ^ (*’> 


‘ ihrf ^ (“> 

*'a 

for mean velocity. Thus, given the surface temperature distribution, 
the total force can be calculated. 


Surface Tension Force 


The total surface tension force tending to close the cavity can 
be calculated as the product of surface tension and the circumfere:.ce 
of the cavity, i.e., 


F = 7T Da 


( 21 ) 


where 

D = cavity diameter 
a = surface tension. 


Hydrostatic Force 

One of the three expected differences in zero-g phenomena as 
compared to one-g is the decrease in hydrostatic pressure exeirted 
by the weight of the molten material on the cavity wall. This is given 
by 


Fj^/A = pgd 


( 22 ) 


where 

d = depth of cavity. 

Thus, the total force is then 

= f PgdD^. 


( 23 ) 
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Electrostatic and Electromagnetic Forces 

The electrostatic force due to a net charge on a mass of material 
(especially low electrical conductors) may be a non -negligible force 
for both e;.periments. The magnitude of this force may be estimated 
for a givgn charge distribution by calculating the electrostatic body 
force (pE) and realizing the tendency to negate surface tension effects. 

The contribution of magnetostriction and electrostriction effects 
(i.e., the variation in electrical properties with changes in density 
due to compression) to the total force may also be estimated if this 
variation in electrical properties (permittivity and permeability) with 
density is known. Also, the Lorentz force (J x 6) is a measure of 
the role of convective charge transport and conduction current of the 
total force. 

bi the absence of one predominant body force (gravity) some 
of tlie mentioned mechanisms may produce a unique effect and the 
relative magnitxides should be established to help clarify the effect 
of the absence of gravity. 


C. ADHESION, COHESION AND SURFACE TENSION EFFECTS 

I lie terms cohesion and surface tension are generally used 
interchangeably while in fact they represent different mechanisms. 

In discussing the differences between the two, it is necessary to 
discuss the concept of intrinsic pressure. Intrinsic pressure is 
generated as a result of the attractive cohesive force producing a 
compression and, since one is a result of the other, they exist at 
tlie same time. It is this intrinsic pressure gradient that defines 
the range of the effects of surface tension itself as surface tension 
is the difference between the intrinsic pressure forces and the co- 
hesive forces. This difference is non-zero only near the surface 
as the coliesive forces are weaker at the surface than they are in the 
bulk and intrinsic pressure is at a maximum. This situation produces 
the largest difference at the surface. The general effect of increasing 
temperature is to superimpose a repulsive force that produces a 
weakened cohesive force and also reduces the intrinsic pressures. 

If the cohesive effects are a stronger function of temperature than 
intrinsic pressure effects, a reduction in the surface tension ac- 
companies an increase in temperature and vice versa. For equal 
dependence, surface tension would remain constant. 


Contact Angle 

The equilibrium conditions of a sensible drop in contact with 
a solid surface have generally been written in terms of surface 
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tension {i.e., the Young-Dupre' form). One derivation of this form 
of the equation can be obtained from a force balance at the interfaces 
(see Figure 2). 



This yields 


SV 


= cose 


+ CT 


LS 


( 24 ) 


A dis Inction must be made, however, between surface tension 
and surface free energy as recognized by Gibbs [l] , restated by 
Shuttleworth [z] and discussed by Pethica and Pethica [3]. Young's 
equation is frequently used to draw conclusions concerning the sur- 
face energy of Solids; but as stated it is objectional from a theoretical 
point of view and refuted by experiment. The main theoretical objec- 
tion lies in the fact that the vertical component (OLysinG) is not 
accovinted for. Some investigators have argued that it is obvious 
that this component is balanced by the elasticity of the solid, by 
gravitation, or by another non-capillary agent. However, if a verti- 
cal plate is envisioned (i.e., unaccounted component now horizontal), 
it is easily seen that this component cannot affect the vertical dis- 
placement of the plate. Also it is known from experiment that poorly 
wetted solid particles float on the surface of a less dense liquid; thus 
the vertical component cannot be compensated for by Ogy - g (see 
Figure 3). These aforementioned sentiments are best described by 
Bikerman [4], as this discussion denies the validity of Young's equa- 
tion from a theoretical and experimental point of view. The work of 
Pethica and Pethica [s] examined the validity of the equation in a 
gravitational field and it was concluded that the equation is invalid. 
Others, as pointed out by Johnson [ 5 ], deny the validity for zero con- 
tact angle. Johnson, however, maintains that Young's equation is 
indeed valid when surface tensions are used instead of surface free 
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FIGURE 3. "FLOATING" SOLID 

energy. He derives relationships based on Gibbsian thermodynamics 
that include the effects of adsorption and gravity. He also discusses 
four major reasons that contribute to these discrepancies, i.e., 

• Assumption that surface tension is numerically 
equal to specific surface free energy 

• Conclusion that a necessary and sufficient condi- 
tion for equilibrium is that the free surface energy 
of the system be a minimum 

• Lack of rigor in defining the differences between 
surface tensions, surface free energies, specific 
surface free energies and surface energies, and 

• Interpretation of a as applied to solids. 

Johnson then clarifies the above misconceptions and derives an equa- 
tion that represents total mechanical equilibrium from which Young's 
surface tension equation results. In support of Johnson's work is 
the work of Li 6 . Here the Young equation is shown to be a terminal 
condition of an isoperimetric problem derived from the principal of 
minimum energy. It is also shown that the contact angle is constant 
for all gravitational field levels. This concept will be employed to 
examine tlie adhesion and cohesion effects in the M551 and M553 
experiments. 
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Spreading and Wetting 


The spreading of a liquid on a substrate has classically been 
defined as the energy released per unit area when the liquid with a 
free surface spreads over the solid. The spreading coefficient is 
given by 


or by 


S = 
c 


'SV "LS 


- a 


LV 


(25) 


(cose - 1) 


(26) 


if Young's equation is substituted. If S < 0 the liqui :issume 

the shape of a spherical segment (no external forces, . i the area 
of contact given by 


A 


1/3 .2. .,2/3 
It ' sin 9 V ‘ 

(1 - COS0) ' ( ^ 


TJl 

) 


(27) 


where 


V = volume of the liquid. 

On the other hand, if > 0, the liquid will spread over the solid 
surface at a rate dependent upon the viscosity of the liquid and the 
surface roughness. 


Work of Adhesion 

Dupre' introduced the concept of reversible work of adhesion 
of liquids to solids and this relation can be written as 

^ ~ ^SV ‘^LV “ ‘^LS 

This equation is simply the thermodynamic expression of the fact that 
the change in free energy of a system is equal to the work needed to 
separate the liquid from the solid. 
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D. LIQUID DYNAMICS 

Formation Time 

The question of the time reeded for rninimum energy formation 
has received considerable attention and is reviewed here to document 
how the calculations have been made. 

Consider a mass of liquid whose motion can be dercribed by the 
full Navier-Stokes equations as given below: 

p = F - VP + p V2 V (d9) 

Performing an order of magnitude analysis with 


V = 

0(r/t) 

D/Dl = 

0(l/t) 

F = 

0(a/r)/: 

V = 

0(l/r) 

v" . 

U(l/r2) 


and keeping track of the resisting and forming forces yields 



pO(r/t^) = ^i^ - 0(P/r).p 0(l/rt) 

i31) 

or 

0(pr/t^) = 0((T/r^) - 0(p/rt) 

(32) 

as 

O (P/r) = O (2o/r2). 

(33) 

Rearranging yields 



0 (ot^) - O (p rt) - O(pr^) = 0 . 

(34) 

Thus 

pr + J fi ^ + 4a p 

'f = 2a 

(35) 


where is the time to formation. For low viscosity liquids such as 
liquid metals, Eq. (35) can be simplified to 
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Damping Time 


After formation, the sphere will continue to oscillate around 
its minimum energv position due to inertial effects. These oscilla- 
tions are reduced as a function of time as a result of internal friction 
due to a non-zero viscosity. 


Consider a system of mass M freely oscillating about its equi- 
librixmi position. Further, it is reasonable to assume a constraining 
force coefficient proportional to the surface tension and a damping 
coefficient proportional to viscosity and radius. This can be expressed 
as 


M^ + C^+Ka = 0 
dt^ 


(38) 


where 

M = mass = (4/3) ir r p 
a = amplitude of response 
t = elapsed time 
C oc rju 
K OC CT . 


The solution to Eq. (38) is 


= Ae 


-tC/2M 


cos ((jt + </>) 


(39) 



J 


J 


1 



where 


A = maximum amplitude 
w = angular velocity 
= phase angle. 

The time required to achieve a certain percent decay can be readily 
found from Eq. (39) as follows: 

2 

t = -2 7 T fn (a/A)oc - ^ /n (a/A) (40) 

O I' 

with the corresponding circular frequency being 



(41) 


E. MELTING AND SOLIDIFICATION TIMES 


Time to Melt Specimen 

From the operating characteristics of a heat source, the time 
needed for melting can be calculated by equating the heat input from 
the source to the change in internal energy plus the radiative loss. 

This yields the following for minimal conduction loss (i.e., M553); 

q = MC^ + Aoe (t'^-T^) (42) 


or rearranging 



(43) 
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Performing this integration yields 


prC 

£n 

. 

• 

1 

(t‘ + Tj^) (t' - T.) 

1/2 

6ac t'3 


(t‘ - Tj^) (T* + T.) 
1 

« 


+ tan (Tj^/T 


) - tan 


(T^T 


where 


I 


.4 


T 



ae 


+ T 


4 

o 


(44) 


(45) 


Upon reaching the melting temperature, the time needed to melt 
can be estimated by equating the available heat (source input less 
radiative loss) to the heat absorbed during formation of a complete 
melt. This gives 


where 


4 - 47t r^ ea = .| rr p . 


M ■ 3 


t,. . = time to melt 
M 

Tj^ = temperature of melting. 


Rearranging yields 


t 


M 


(4/3) 7T Aj r^ p 

” 2 
o - 4 )t r € a T., 
M 


(46) 


(47) 


Thus, the total time for heating and melting is 

*^total ^ ^ 

In the derivations for both of these times, it was assumed that 
isothermal conditions prevailed; i.e., the sphere is at uniform tem- 
perature. This would be realistic for materials of relatively high 
thermal conductivity. However, a "warm up" time is associated 
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with low k materials as a result o£ finite heat conduction and this 
time can be estimated by an order of magnitude analysis of Fourier's 
heat conduction equation, 

v2T = .i|I (49) 

with 

= 0(l/r)2 (50) 

a = thermal diffusivity = k/pC 

t = conduction time 
c 

Thus 

0(T/r^) = ^0(T/t^) (51) 

or rearranging 

0(l/r^) = 0(l/at^). (52) 

This relation shows that the thermal diffusivity is a measure of the 
time required to heat a material to a given temperature level and is 
seen to be proportional to the square of the conducting path and 
reciprocal transfer time. A final result can be obtained from Eq. 

(52) by solving for this time, 

t^ = rV« = p Cp/k (53) 

In a transient heating process the thermal capacity of the 
material governs the amount of energy absorbed and the thermal 
conductivity dictates the rate. Thus, for large k the time to reach 
the melting temperature is governed by Eq. (44) while for small k, 

Eq. (53) applies. A maximum time (upper bound) can be estimated 
by the addition of all three times (Eqs. (44), (47) and (53)); however, 
in practice a tradeoff should be made between Eq. (44) and Eq. (53), 
depending on the value of the conductivity and the magnitude of con- 
vection present. 

The maximum iemperature capable of a given heat source can 
be readily obtained from a balance of heat input from source and the 
loss due to radiation. This can be expressed as 

q = 4^” r^ € CT T"^ (54) 
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or 


T 

max 



(55) 


Solidification Time 


The time needed for a liquid sphere to solidify can be estimated 
by equating the heat of fusion to the time integral of the radiated dis> 
sipation rate; i.e., 


M 



Ae a (T^ 


'4 

- T dt 


(56) 


where 

= heat of fusion 
A = surface area 
€ = emissivity 

a = Stefan-Boltzmann constant 
T = solidification temperature 

I 

T = surrounding temperature. 

'4 4 

Since T « T , Eq. (56) can be reduced to 


4 

3e oT^ 


(57) 
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SECTION III. M551 METALS MELTING EXPERIMENT 
A. GENERAL 

During the course of the investigation all process phenorrena 
have been considered with particular attention given to adhesion - 
cohesion studies with emphasis on the fluid dynamics of the molten 
metal. This was accomplished by establishing all expected varia- 
tions in terrestrial versus space processing occurring as a result 
of 

• The dominance of capillary motion (instead of gravity- 
induced flow) 

• Minimal segregation (sedimentation), and 

• A decrease in hydrostatic pressure gradient. 

These mechanisms are all hydrodynamical in nature; thus gravity 
has no direct affect on grain structure or any other properties of 
the solidified material. The gravity effects occur as a result of 
differing fluid motion . 

The following items are considered the most significant contri- 
bution to the M551 experiment analysis: 

• Establishment of the role of vaporization in the 
cavity shape 

• Performance of three-dimensional temperature 
calculations for weld and dwell modes of melting 
for multi -component variable property material, 
and 

• Establishment of magnitude and direction of gravity 
during actual Skylab mission to allow meaningful 
conclusions to be drawn from mic’ ographs about 
zero-gravity effects. 

Analysis of the aforementioned results has produced the fol- 
lowing projections of effects of zero-gravity. 

• Since all fluid mechanic phenomena require a finite 
time to occur (due to inertia considerations) an 
effect due to gravity during the weld portion of the 
experiment may not occur (except for a more spheri- 
cal bead when beading occurs) as large cooling rates 
are available via conduction. 
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• A more homogeneous dwell structure will occur 
due to lack of sedimentation. 

• The liquid cavity shape will be determined by 
vaporization effects and should be nearly symmetric 
(different temperature gradients do exist in positive 
and negative radial directions). 

• If degassing is the cause of the splattering as noted 
in color films of the process, porosity may result 
in the zero-gravity dwell as the driving force for 
the pockets to surface is less. 


B. THERMAL HISTORY RESULTS 

Thermal histories corresponding to the experimental test con- 
ducted at NASA -MSEC with the M512 Facility have been calculated. 
The stainless steel disk was discretely divided into 289 nodal points 
shown schematically in Figure 4. The operating procedure consisted 
of the weld beginning at node 4 and continuing arovmd the weld circle 
through all four quadrants. The weld portion ends at a time of ap- 
proxim;.tely 19 seconds with the dwell beginning on node 316 at 19.5 
seconds and continuing xintil 34.5 seconds (15-second duration). The 
input power was 1 kW with a 75% conversion efficiency assumed for 
the calculations. Table 1 defines the eb location as a function of 
time for both the weld and dwell modes of operation. 

Figure 5 shows the nodes on the weld circle for the entire 
quadrant of the thin section (0.02-inch thickness). This result indi- 
cates that fxill penetration will occur with an average of 500°C super- 
heat. Notice also as the dwell is begun at 20 seconds, no effect is 
seen at this location. Figure 6 shows the radial distribution of tem- 
perature corresponding to the center of the quadrant. Note that the 
nodes surrounding the weld circle do not melt as the maximum tem- 
perature is 800°C. Also a nearly symmetric radial temperature 
gradient about the weld circle is seen to exist. Figure 7 gives the 
weld circle distribution of temperature through the entire quadrant 
of 0.05-inch thickness. Again, full penetration is predicted with 
300°C of superheat. Figure 8 reveals the same results as those 
for the thin section, namely, no melting except on the weld circle 
and again near symmetry existing. 

The situation in the third quadrant consisting of a 0.125-inch 
thickness section is shown in Figure 9- For this beam power (1 kW) 
these results indicate that melting is confined to near the surface. 

In fact the average temperature over one third of the thickness is 
below the melt temperature indicating that less than one-third pene- 
tration occurs. Figure 10 shows the profiles just below the surface 
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Table 1 

ELECTRON BEAM POSITION VERSUS TIME iDISK ROTATION) 


Node 

Time 

(sec) 

Node 

Time 

(sec) 

Node 

Time 

(sec) 

Node 

Time 

(sec) 

4 

0 

■1 

5.5805 

204 

11.161 

■1 


10 

0.62005 

wm 

6.2005 

204 

11.657 

mm 

18.454 

16 

1.2401 


6.8206 

210 

12.773 

316 

19.570 

22 

1.8602 

wm 

7.4406 

216 

13.889 

316 

34.570 

28 

2.4802 

■a 

8.0607 

222 

15.005 



34 

3.1003 

134 

8.6807 

228 

16.121 



40 

3.7203 

140 

9.3008 





46 

4.3404 

146 

9.9208 





52 

4.9604 

152 

10.541 
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FIGURE 4. NODAL POINT ARRANGEMENT FOR THERMAL CALCULATIONS 
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FIGURP' 5. THERMAL HISTORY FOR 0.02-INCH STAINLESS STEEL DISK 
(NODES 4, 10, 16, 28,40 AND 52) 
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FIGURE 6. THERMAL HISTORY FOR 0.02-INCH STAINLESS STEEL 
DISK (NODES 3’. THROUGH 36) 
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FIGURE 7. THERMAL HISTORY FOR 0.05-INCH STAINLESS STEEL 
DISK (NODES 104, 110, 1 16, 128, 140 AND 152) 
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FIGURE 8. 'iHERMAT. HISTORY FOR 0,05-INCH STAINLESS STEEL 
DISK (NODES 131 THROUGH 136) 
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FIGURE 10. THERMAL HISTORY FOR 0.125-INCH STAINLESS STEEL 
DISK (NODES 249 THROUGH 254) 
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node as shown in the previous figure. Exactly the same trend is 
observed with lower maximum temperatures. The radial profile 
for the bottom layer is shown in Figure 11, indicating precisely the 
same trend as before. The thickest quadrant results are shown in 
Figure 12 for the upper surface where again melting is seen to occur 
only near the surface. Here, an asymmetric situation develops about 
the weld circle with a substantial difference in temperature in the 
positive and negative radial direction?. As the center of this quadrant 
is reached the dwell portion of the cycle is initiated as depicted in 
Figures 13 and 14. It is seen that nearly full penetration is predicted 
(no allowance for eb cavity) and a substantial melt region is outlined. 
Witltin a few sec onds of eb cutoff, nearly isothermal conditions result 
near the dwell region as indicated by the temperatures at 40 seconds. 

A comparison between these calculated results and those ob- 
tained by experiment are shown in Table 2 where the maximum tem- 
peratures are given in each case. The calculated results in the weld 
case are consistently lower than the corresponding experimental 
values. As previously stated, a 1 kW beam power (50 mA and 20 kV) 
was assumed with an efficiency of 75% giving rise to an input power of 
0.75 kW. If the actual input power was slightly higher, this could 
explain the trend. Also some of the maximum temperatures for the. 
dwell occur after 40 seconds. This accounts for the dwell tempera- 
tures iii the table being listed as greater than some level (i.e., > 50). 
This means that the maximum temperature will occur later in time 
than the cutoff of the plot. In general, however, agreement is fairly 
good. Temperature histories were also calculated for aluminum and 
tantalum disks, and the results are given in Appendix B. 


C. BEADINC. CONSIDER ATIONS 

Upon examination of the M5 51 disks, it has been noted that a 
beading occurs behind the electron beam as the beam moves through 
the molten material. At present, no detailed analysis has been con- 
ducted to explain this phenomenon; however, some consideration has 
been given to understanding the mechanism. In order to present a 
hypothesis explaining the proposed mechanism some backgrotmd 
information is presented. 

The movement of the electron beam cavity through the molten 
metal resembles the motion of an infinite cylinder (especially if 
full penetration occurs) through a fluid. It has long been known that 
such a body leaves in its wake a regular pattern of vortices which 
move alternately clockwise and counterclockwise and is known as a 
Karman vortex street. A schematic is presented in Figure 15. 
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FIGURE 11. THERMAL HISTORY FOR 0.125-INCH STAINLESS STEEL 
DISK (NODES 279 THROUGH 284) 
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FIGURE 12. THERlvlAL HISTORY FOR 0.25-INCH STAINLESS STEEL 
DISK {NODES 301 THROUGH 306) 
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figure 14. THERMAL HISTORY FOR 0.25-INCH STAINLESS STEEL 
DISK (NODES 310, 322, 328, 316, 346 AND 376) 
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FIGURE 15. KARMAN VORTEX STREET (DIAGRAMMATIC); 
STREAMLINES DRAWN IN A SYSTEM OF 
COORDINATES MOVING WITH THE VORTEX 
STREET 


There exists a distinct frequency at which the vortices on one 
side of the wake are shed that depends only on the Reynolds number 
of the fluid motion and is given by the Strouhal number, i.e., 



where f :s the frequency, D is the diameter of the cavity and v is 
the veloc ity of the beam in the liquid metal. 

From experiment, the relationship of the Strouhal number to 
tlie Reynolds number has been determined and is shown in Figure 
16. i his shows that a critical Reynolds number exists below which 
no shedding has been observed. 

Using these concepts the theoretical number of beads per unit 
length lias been determined for the M551 materials. It should be 
mentioned that no attempt has been made to account for the solid 
material (edges) in the cut or the extent of melt in the radial direc- 
tion of the disks (which is a function of the thermal conductivity). 
These results are shown in Table 3. The surface tension is included 
as the magnitude indicates the tendency of the vortices to bead before 
solidification occurs. 
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FIGURE 16. STROUHAL NUMBER VS REYNOLDS 
NUMBER 


TABLE 3 

BEAD POPULATION 


M551 

Materials 

Reynolds 

No. 

Strouhal 

No. 

Shedding 

Freq. 

(beads/ 

sec) 

Beads 
per cm 
(on each 
side) 

Total 

Beads 

(both 

sides) 

Surface 

Tension 

(dynes/ 

cm) 

321 Stainless 
Steel 

55 

0.125 

0.8 

m 

1.0 

1750 

2219 Alumi- 
num 

86 

0.155 

1.0 

0.6 

1.2 

737 

Tantalum 

207 

0.19 

1.2 

0.75 

1.5 

2150 


Some qxialitative observation can be made concerning these 
results. First both aluminum and stainless possess Reynolds nvim- 
bers near the critical values; thus it would not be surprising if some 
of the effects not accounted for could have a significant effect on the 
population of the beads if not controlling whether or not beading will 
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I occur. Tantalum, however, possesses the largest Re/nolds number 

and value of surface tens!' t thus it is expected that these properties 
produce the greatest probuoility of beading and will give rise to the 
greatest population of beads. 

! 

I D. FLIGHT RESULTS 

f Shown in Figures 17 through 22 are photographs taken of stain* 

r less steel, aluminum and tantalum metal disks alter being subjected 

to electron beam iinpingement in the M551 metals melting experi- 
ment. In these photographs, Skylab flight test specimens are com- 
pared to ground-based test specimens. The most striking phenome- 
non appearing in these photographs is the "beading" effect which 
occurs simultaneously with total penetration in the stainless steel 
and tantalum specimens. It was speculated in Section III-C that this 
beading effect may be related to vortex shedding. Theoretical bead 
spacing along the electron beam track was predicted based on cal- 
culated vortex shedding frequencies. Average bead spacing distances 
were estimated for both flight test and ground based test specimens 
by inspection of the photographs in Figures 17 through 22. These 
experimental values are compared to the predicted values (Table 3) 
in Table 4. 


TABLE 4 
BEAD SPACING 


M551 

Materials 

Bead Spacing, cm 

Theoretical 

Flight 

Test 

Ground Based 
Test 

321 Stainless Steel 

1.0 

0.356 

0.418 

22 1 9 Aluminum 

0.8 

— 

— 

Tantalum 

0.7 

1.57 

0.94 


The experimental bead spacing on the stainless steel and tantalum 
specimens is seen to be in good qualitative agreement with the pre- 
dicted value for both flight and ground based tests. No beading is 
exhibited on the aluminum specimen, however, and the spacing on 
the tantalum specimen is larger than the predicted value. Consider- 
ation was given to the possibility that the departure from vortex 
sliedding theory might be due to "waill" effects resulting from the 
limitation on the extent of melting. The presence of confining parallel 
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FIGURE 18. STAINLESS STEEL DISK AFTER M551 GROUND-BASED TESTS 
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FIGURE 19. ALUMINUM DISK AFTER M551 FLIGHT TESTS 
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FIGURE 22. TANTALUM- DISK AFTER M551 GROUND-BASED TESTS 
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walls will tend to increase the drag force on the cylinder due to the 
increased fluid velocity in the space between cylinder and wall. This 
may be considered equivalent to an increased "effective" Reynolds 
number for drag on the cylinder. Hypothesizing that the Strouhal 
nvimber relates to this effective Reynolds number according to the 
relationship in Figure 16, the vortex shedding frequency should be 
increased by the presence of the confining walls. The results in 
Table 4, however, indicate the opposite trend. Wall effects, there- 
fore, do not accovint for the observed departure from vortex shedding 
theory. Also, nearly identical bead spacing is noted for both flight 
and ground-based tests with the stainless steel specimens, while 
spacing on the tantalum specimens is larger for flight than for ground- 
baced test. No definite trend can be established, therefore, in the 
effect of reduced gravity on bead spacing. 

An indication of total penetration depth can also be obtained by 
inspection of the photographs. The plate thickness is smallest at the 
starting point of the electron beam track and increases with angular 
displacement along the track. Total penetration is seen to occur on 
all plates beginning at the beam starting point and ceasing at angular 
positions where the plate thickness becomes greater than the penetra- 
tion depth. The penetration depth, therefore, should be approximately 
equal to the plate thickness where total penetration ceases. The angu- 
lar positions where total penetration ceases and the corresponding 
plate thicknesses are given in Table 5 along with theoretical penetra- 
tion depths determined according to a method presented in [t^ . The 


TABLE 5 

OBSERVED AND THEORETICAL 
PENETRATION DEPTHS 



Flight Test 

Ground Based 
Test 


M551 

Materials 

Plate 

Position 

(deg) 

Plate 

Thickness 

(cm) 

Plate 

Position 

(deg) 

Plate 

Thickness 

(cm) 

Theoretical 

Depth 

321 Stainless 
Steel 

75 

0.104 

70 

0.915 

0.195 

2219 Alumi- 
num 

115 

0.211 

108 

0.193 

0.323 

Tantalum 

98 


125 

0.762 

0.120 


f JfTiWf? mob BIOTS NO? mm 


1 













( theoretical method is based on an energy balance in which the beam 

I power is equated to energy absorbed in melting and lost to thermal 

I conduction and re-radiation. The theoretical depths are seen to be 

> within about a factor of two of the experimenttil results in all cases. 

! No definite trends are noted in the differenceti between flight test 

I and ground-based test results. 
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SECTION rv. M553 SPHERE FORMING EXPERIMENT 


A. GENERAL 

During the course of the investigation all process phenomena 
have been considered with particular attention given to adhesion- 
cohesion studies with emphasis on the fluid dynamics of the molten 
metal. This was accomplished by establishing all expected variations 
in terrestrial versus space processing occurring as a result of; 

• The dominance of capillary motion (instead of gravity- 
induced flow) 

• Minimal segregation (sedimentation), and 

• Decrease in hydrostatic pressure gradient. 

These mechanisms are all hydrodynamic in nature, thus gravity has 
no direct effect on grain structure or any other property of the 
solidified material; the gravity effects occur as a result of differing 
fluid motion . 

The following items are considered the most significant contri- 
butions to the M553 experiment analysis: 

• Identification of physical forces affecting the tra- 
jectory of the molten sphere 

• Development of unifying calculational technique to 
predict trajectory of molten spheres during free- 
float condition 

• Development of three-dimensional thermal model 
to predict temperature history of sphere during 
melting 

• Establishment of the role of vaporization on 
temperature history and motion of sphere 

• Establishment of the magnitude and qualitative 
motion of surface dynamics, and 

• Performance of detailed analysis of KC-135 flight 
test data and ground-based experimental tests. 

Analysis of the aforementioned results has produced the fol- 
lowing projections of effects of zero-gravity: 
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• Less heterogeneous nucleation will occur due to 
foreign particles as segregation effects will be 
minimal. 

• If the material has any gas content that comes out 
of solution during solidification, an increase in 
porosity may result as there will be minimal 
tendency for the gas pockets to surface (except for 
convection effects). 

• For the materials with high vaporization rates 
(Ni-1 2% Sn, and Ni*30% Cu), a substantial increase 
in surface cooling occurs — thus greater nucleation 
and surface tension gradients. 

• Even though the cooling rate is dependent on the 
amount of convection, the cooling rate will be deter- 
mined by radiation considerations only as the sphere 
reaches isothermal superheat conditions in approxi- 
mately a second. This in turn dictates the solidifica- 
tion rate. 

• Since solid formations in polycrystalline materials is 
highly dependent on the distribution of foreign particles 
and the liquid motion, the associated differences should 
be discernible (homogeneous nucleation). 

As discussed previously, the only expected differences in zero- 
g versus one-g processing are the dominance of capillary motion 
instead of gravity-induced flow, the absence of segregation and a 
decrease in hydrostatic pressure. There are, however, many in- 
direct effects (advantages) due to a lack of gravity. In general, the 
factors affecting the final sample properties, composition and struc- 
ture can be divided into three major categories, i.e., 

• Nucleation 

• Growth, and 

• Reaction kinetics. 

The absence of gravity might allow better control of the independent 
variables affecting each category; this occurs directly through the 
three aforementioned fluid dynamic effects and indirectly through 
other mechanisms that are a function of fluid motion (gas content, 
composition, foreign particles, imperfections, cooling rate, nuclea- 
tion rate, etc.). 
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B. M512 VACUUM CHAMBER GEOMETRY 

The geometry of the M512 chamber is shown in Figure 23 along 
with the coordinate system chosen. The pinwheel containing the speci- 
men to be melted is mounted in the chamber in the y-z plane with the 
center of the specimen to be melted positioned at coordinates (-0.292, 
-0.750, -0.433) as shown in Figure 24. Figure 24 also defines the 
camera locations used in the KC-135 aircraft and ground-based 
(camera 1) tests with the associated line-of-sight geometry to account 
for mirror effects. The corresponding film frame projections are 
depicted in Figure 25 in the local 4* C coordinate system. 


C. M553 THERMAL CONSIDERATIONS 

The theoretical foundation for computing the three-dimensional 
temperature history for a multi- component system has been pre- 
sented in Section II. This subsection presents the results of calcula- 
tions for the M553 materials of interest. First consider the time 
required for the material to become superheated and cool to the melt 
temperature for the operating procedures used in the actual Skylab 
mission. This time can be determined from the temperature history 
calculated by the Lockheed Thermal Analyzer which includes allow- 
ances for; 

• A three-dimensional network 

• Conduction heat loss 

• Radiation heat loss 

• Vaporization heat loss 

• Variable properties (with temperature and phase), and 

• Variable heat source. 

The spheres are discretely modeled by 69 nodal points with each 
of the corresponding nodal temperature histories computed. Allow- 
ances are included for sting melting and automatic electron beam 
cutoff. The nodal arrangement is shown for the inner layer in Figure 
26, with the outer layer and electron beam impingement location 
given in Figure 27. 

The thermal history for pure nickel is shown in Figures 28 
through 31. It is seen that in less than one second after eb cutoff 
the sphere is nearly isothermal in each layer. The time correspond- 
ing to the onset of solidification is shown in Figure 28 and is approxi- 
mately 10.5 seconds. This time must be added to the time needed for 
solidification to determine the necessary free-float time for con- 
tainerless solidification. The corresponding temperature histories 
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FIGURE 24. RELATIVE GEOMETRY OF M512 CHAMBER AND CAMERA 



Camera 1 Film Frame 



FIGURE 25. PROJECTION OF x'y'z' AXES ONTO FILM FRAME COORDINATE 
SYSTEM (POSITIVE NORMAL, IS INTO PLANE) 
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FIGURE 27. , THE RELATIVE POSITION OF ELECTRON BEAM TO SPHERE 
AND 1 HFPMAL MODEL FOR OUTER LAVER 
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FIGURE 28. NICKEL TEMPERATURE HISTORY (NODES 1,2,3,4AND8) 
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FIGTIRE 30. NICKEL TEMPERATURE HISTORY (NODES 36, 37. 38, 42 AND 46) 
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FIGURE 31. NICKEL TEMPERATURE HISTORY (NODES 50, 54, 58, 62 AND 64) 
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for Ni-12% Sn, Ni-1% Ag and Ni-30% Cu are shown in Figures 32 
through 35, 36 through 39 and 40 through 43, respectively. The 
essential features of these results are as follows: 

• Vaporization forces can be computed 

• Time above melt temperature can be determined, 
and 

• Superheat can be determined. 

These items either affect the free-float time directly or the time 
duration needed for solidification. The solidification time including 
both vaporization and radiation losses for all four materials has been 
computed based on the methodology of Section II. These resvilts are 
summarized in Table 6 for the M553 materials. These calculations 
indicate that from 32 to 48 seconds are necessary for complete 
solidification to occur. 

Unfortunately, provision was not made during the Skylab flight 
tests for temperature measurements. Grovind test measurements 
were maae, however, and reported in Reference . The measured 
temperaa’res represent a spatial average over the observed sphere 
surface (near the impingement point) and a time average for a period 
on the order of one or two seconds after electron beam cutoff. For 
the pure nickel specimens, the maximum measured temperature was 
about 1500°C for the sting-type specimens and about 1450°C for the 
release specimens. Figures 2C and 29 show predicted "average" 
temperatures of about 1600°C for the one or two second time period 
after beam cutoff. An emittance correction for the sample could have 
raised the measured temperature by about 100°C thus giving good agree- 
ment between the predicted and measured results [8]. The maximum 
measured temperature for the nickel- 12% tin specimens was about 
1400 C, which after emittance corrections could be increased to about 
1500*^C. This is about 100°C less than the 1600°C predicted average 
indicated in Figures 32 through 36. The maximum measured tempera- 
ture for the nickel- 1% silver specimens was also about 1400°C with 
possible emittance correction to 1500°C, and this also compares to a 
predicted value of about 1600°C indicated in Figures 36 through 39. 

For the nickel-30% copper specimens, a maximum measured value of 
about 1350°C compares to a predicted value (from Figures 40 through 
43) of about 1425^G, thus giving good agreement if emittance corrections 
increased the measured values by 50 to 100°C. 
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FIGURE 32. NICKEL-12% TIN TEMPERATURE HISTORY (NODES 1,2, 
3. 4 AND 8) 
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FIGURE 33. NICKEL- 12% TIN TEMPERATURE HISTORY (NODES 12, 18, 
20, 24 AND 32) 
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FIGURE 34. NICKEL- 12% TIN TEMPERATURE HISTORY (NODES 36, 
37, 38, 42 AND 46) 
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FIGURE 35. NICKEL-12% TIN TEMPERATURE HISTORY (NODES 50, 
54, 58, 62 AND 64) 
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FIGURE 36. NICKEL- 1% SILVER TEMPERATURE HISTORIES (NODES 
1,2, 3,4 AND 8) 
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FIGURE 37. NICKEL. 1% SILVER TEMPERATURE HISTORIES (NODES 
36, 37, 38, 42 AND 46) 
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FIGURE 33. NICKEL- 1% SILVER TEMPERATURE HISTC RIES (NODES 
36, 37, 38, 42 AND 46) 


3 



TEMPERATURE CDEG C 



74 


TEMPERATURE CDE6 C 


1 



5 


TEMPERATURE (DEG C. 





TEMPERATURE CDEG C 


■■■r/ABWHail 


NODE 

36 

37 
- 38 

42 

46 


SYNB 

□ 

A 

O 

0 

a 


am 


mmm 

mam 

BBarnaurnaanm 


TIME CSEC) 


FIGURE 42. NlCKEL-30% COPPER TEMPERATURE HISTORIES (NODES 
36, 37. 38, 42 AND 46) 
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FIGURE 43. NICKEL,-30% COPPER TEMPERATURE HISTORIES (NODES 
50, 54, 58, 62 AND 64) 
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D. M553 TRAJECTORY ANALYSIS 


One of the most important aspects of the M553 experiment is 
the motion of the sphere after deployment since this determines what 
free -float time will exist. In order to predict the motion, a trajectory 
program has been developed that includes allowances for: 

• Electron beam force 

• Deployment velocity (spring) 

• Vrporization force based on 3-D temperature 
i .story 

• S: /^lab orbit considerations 

• M512 position considerations, and 

• Allowances for additional forces as independent 
subroutines (KC-135 g-level data). 

In general a solution to the equation of motion 



(58) 


is obtained numerically via Simpson's rule subject to the following 
boundary conditions: 


R - 0 and 


dR 

dt 


= V at 
o 


t = 0 


( 59 ) 


Integration of the equation of motion twice yields 


R 


t rj 



o o 


(60) 


which produces the desired trajectory b' chamber coordinates centered 
on the specimen. This computed trajectory is then transformed to 
camerr coordinates to shew the path as seen in the film, i.e.. 


R'(t) = [A]S(t) (61) 

where R (t) is the projected trajectory and [a] is the transformation 
matrix obtained from the hardware. P or camera 1 
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and for ramera E 



.932 

-.203 

.301 


.661 

-.224 

-.716 


-.363 

-.521 

.773 


-.183 

-.974 

.136 


0 

-.829 

-.560 


-.728 

.041 

.685 


( 62 ) 


(63) 


To demonstrate the concept, consider the accelerometer c?->*^a 
in Table 7. These data were obtained by orthogonal accelerometers 
placed on board the July 1972 KC-135 research aircraft flight. The 
resulting directions of the residual g-level are shown in Figure 44 
transformed uito both camera planes. 

The tabular gravity levels appear in Table 8 in film plane 
( oordinates. The resulting float time as a function of the magnitude 
of tlie gi-avity level is shown in Figure 45 for a spherical outer 
boundary 20 cm from specimen release. It is seen that float times 
of 1 to 3 seconds can be expected in the KC-135 flights fcr g-levels 
of 6(10)-2 to 0.6(1P)"^. 

The eb force was discussed in Section II and Table 9 gives the 
magnitude of this force for two eb currents. To establish the motion 
as a result of the eb force, the conversion efficiency and eb cutoff 
time must be known. This time corresponds to the time delay asso- 
ciated with automatic shutdown of the eb due to sting melting. 

Original estimates were from 10 to 20 milliseconds; however, high 
speed movies (1000 frame/sec) take., during the M553 ground-based 
test at MSFC on 20 Octooer 1972 show conclusively that longer tin;es 
exist. During their test the spring loads were varied on the shutdown 
mechanism to experimentally establish the time interval the eb is 
on after melting occurs. It was found that for a 150-g spring load a 
cutoff time of 410 milliseconds was measured while 220 milliseconds ' 
were observed for a 200- g loac*. This cutoff time translates into a 
free-float time of approximately 12 seconds for a 150-g load and 23 
seconds for the 200-g load. Thus the expected free-float time is 
drastically reduced from the 250 seconds corresponding to the 
original 20 milliseconds estimate. 

An analysis was performed to determine if the cohesive energy- 
generated kinetic energy is sufficient to cause deployment in the 
absence of external forces. This was done to allow initial conditions 
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Table 7 

ACCELEROMETER DATA FROM JULY KC-135 FLIGHT 


Index 

Specimen 

(g) 

a (10)^ 
V 

(g) 

n (10)^ 

(g) 

lg|(10)^ 

(g) 

1 

10 

1.0 

-1.0 

3.0 

3.4 

2 

9 

0.6 

0.3 

0.3 

0.65 

3 

8 

-0.4 

0.3 

4.0 

4.05 

4 

14 

0.0 

0.2 

2.0 

2.02 

5 

11 

-1.0 

-0.7 

2.0 

2.34 

6 

7 

0.6 

-0.4 

1.0 

1.19 

7 

6 

0.6 

-0.1 

1.5 

1.65 

8 

5 

1.4 

-1.0 


4.82 

9 

2 

1.6 

-2.2 


3.14 

10 

15 

-0.4 

-0.4 

-l.O 

1.15 

11 

13 

0.2 

-0.4 

5.8 

5.82 
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Table 8 

TRANSFORMED KC-135 ACCELEROMETER DATA 


Camera 1 

Index 

Specimen 


^1 


1 

10 

1.290 

-2.170 

-2.150 

2 

9 

0.357 

-0.506 

0.214 

3 

8 

-0.482 

-3.390 

-2.125 

4 

14 

-0.073 

-1.762 

-0.964 

5 

11 

-0.678 

-1.090 

-1.960 

6 

7 

0.61 1 

-0.722 

-0.718 

7 

6 

0.596 

-1.313 

-0.736 

8 

5 

1.670 

-3.490 

-2.868 

9 

2 

2.195 

2.666 

-0.019 


15 

-0.228 

1.119 

0.130 

1 1 

13 

0.331 

-4.641 

-3.492 


Camera 2 

1 

Index 

Specimen 

^2 

1 

^2 

^2 

1 

10 

-1.339 

0.874 

-2.906 

2 

9 

0.057 

-0.392 

-0.523 

3 

8 

-3.230 

-0.038 

-2. 411 

4 

14 

-1.492 

-0.112 

-1.342 

5 

11 

-1.989 

0.988 

-0.748 

6 

7 

-0.324 

0.319 

-1.097 

7 

6 

-0.677 

0.025 

-1.470 

8 

5 

-2.167 

0.846 

1.716 

9 

2 

2.99 

1.716 

0.134 


15 

0.536 

0.438 

0'.917 

11 

13 

-4.016 

0.584 

-4.168 
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Table 9 

TOTAL FORCE DUE TO ELECTRON BEAM 

I 2 m V 

_ t ^ e 


I 

e 

m 

e 

V 

1 

F 

e 

F 

e 

(amp) 

(coulomb) 

(kg) 

(V) 

(N) 

(dyne) 

0.05 

1.6(10)‘^'^ 

-30 

0.9(10) 

2(10)^ 

2.38(10)“^ 

2.38 

0.10 

1 1 

1.6(10)'^^ 

-30 

0.9(10) 

2(10)^ 

4.76 (10)"^ 

4.76 
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FIGURE 45. APPROXIMATE FLOAT TIME AS A FUNCTION OF GRAVITY 
LEVEL 
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to be calculated. To reduce the complexity of the problem, the fol- 
lowing assumptions were made: 

• Zero gravity field 

• No sting melting 

• Sting pulled down before breaking 

• Perfectly spherical free surface, and 

• Contact angle of 120 degrees. 


An energy balance was made between the initial and .‘'inal states 
and from this the kinetic energy is found as 


KE 


LV 



+ 


This kinetic energy was calculated as a function of penetration 
(H) of liquid into orifice. The value of H simulates the effect of the 
spring in pulling the liquid down via adhesion between liquid and solid 
sting. The results of this calculation are 


Penetration, H 

Kinetic Energy 

(cm) 

(ergs) 

0.0 

0.0 

0.1 

62.7 

0.5 

293.0 

l.O 

597.0 


The value of the kinetic energy must then be greater than the 
work of adhesion (or energy of adhesion) for deployment to occur. • 
The work of adhesion was found to be 265 ergs; thus for depressions 
cf more the 0.5 cm deployment will occur due to cohesive effects 
alone. These criteria are included in the trajectory computer pro- 
gram to allow initial conditions for the numerical integration to be 
determined. 

Tables 10 and 11 show the surface temperature distribution 
at 3.0 and 3.5 seconds, respectively, for Ni and Ni-12% Sn. These 
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Table 10 

TEMPERATURE DISTRIBUTION AT 3.0 SECONDS 



Nickel 

Nickel - 12% Tin 

Node 

Temp. 

Vapor Pressure 

Temp. 

Vapor Pressure 

No. 

(‘'C) 

(mm Hg/cm^) 

(°C) 

(mm Hg/ cm2) 


2036 

.3959 

2013 

.7001 


1463 

.0006 

1537 

.0069 

38 

2193 

.464 

2157 

.6965 

39 

1900 

.0395 

1903 

.0962 

40 

1728 

.0067 

1757 

.0253 

41 

1677 

.0037 

1716 

.0168 

42 

1677 

.0037 

1716 

.0168 

43 

1728 

.0067 


. .0254 

44 

1900 

.0400 

1903 

.0963 

45 

2193 

.4644 

2157 

.6965 

46 

1801 

.0202 

1815 

.0602 

47 

1695 

.0063 

1731 

.0266 

48 

1609 

.0022 

1663 

.0132 

49 

1573 

.0014 

1634 

.0096 

50 

1573 

.0014 

1634 

.0096 

51 

1609 

.0022 

1663 

.0132 

52 

1695 

.0063 

1731 

.0266 

53 

1801 

.0201 

1815 

.0602 

54 

1633 

.0030 

1680 

.0159 

55 

1584 

.0016 

1642 

.0105 

56 

1537 

.0009 

1605 

.0069 

57 

1513 

.0006 

1586 

.0056 

58 

1513 

.0006 

1586 

.0056 

59 

1537 

.0009 

1605 

.0069 

60 

1584 

.0016 

1642 

. 0105 

61 

1633 

.0020 

1680 

.0158 

62 

1535 

.0006 

1602 

.0049 

63 

1515 

.0005 

1585 

.0041 

64 

1492 

.0003 

1567 

.0033 

6 5 

1480 

.0003 

1557 

.0030 

66 

1480 

.0003 

155'’ 

.0030 

67 

1492 

.0003 

1567 

.0033 

68 

1515 

.0005 

1585 

.0041 

69 

1535 

.0006 

1602 

.0049 
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Table 11 


TEMPERATURE DISTRIBUTION AT 3.? SECONDS 


Node 

No. 

Nickel 

Nickel 

-12% Tin 

T emp. 
('^C) 

Vapor Pressure 
(mm Hg/cm^) 

Temp. 

(®C) 

Vapor Pressure 
(mm Jig/cm^) 

36 

1689 

.01249 

1678 

.0332 

37 

160, 

.00429 

l624 

.0185 

38 

1685 

.0041 

1676 

.0111 

39 

1679 

,0038 

1674 

.0109 

40 

J 6 1 2 

.0035 

1671 

.0105 

41 

.667 

.0033 

1670 

.0103 

42 

1667 

.003 3 

1670 

.0103 

43 

1672 

.0035 

1671 

.0105 

44 

167? 

.0038 

1674 

.0109 

45 

1685 

.0041 

1676 

.011 1 

46 

1675 

.0050 

1672 

.0144 

47 

1668 

.0046 

1 669 

.0140 

48 

1658 

.0041 

1665 

.0135 

49 

1652 

.0038 

1662 

.0131 

50 

1652 

.0038 

1662 

.0131 

51 

1658 

.0041 

1665 

.0135 

52 

1668 

.0046 

1669 

.0140 

53 

1675 

.00 50 

1672 

.0144 

54 

1657 

.0040 

1663 

.01 J5l 

55 

1651 

.0037 

1660 

.0128 

56 

1642 

.0033 

1657 

.0123 

57 

1636 

.0031 

1654 

.0200 

58 

1636 

.0031 

1654 

.0200 

59 

1642 

.0033 

1657 

.012 ^ 

60 

1651 

.0037 

1660 

.OliiS 

61 

1657 

.0040 

1 663 

.0131 

62 

1634 

.0022 

1649 

.0083 

63 

1630 

.0021 

1648 

.0082 

64 

1624 

.0020 

1645 

.0080 

65 

1621 

.0019 

1644 

.0078 

66 

1621 

.0019 

1644 

.0078 

67 

1624 

.0020 

1645 

.0080 

68 

1630 

.0021 

1648 

.0082 

69 

1634 

.0022 

1649 

.0083 
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data were used in conjunction with the Langmuir theory discussed in 
Section II to calculate the net vaporization force as a function of time. 
A sample calculation is presented in Table 12 for both Ni and Ni-12% 
Sn for the first se 'ond after eb cutoff. It is seen that between 3.5 
and 4.25 seconds the vaporization force is of the same order of magni* 
tude as the eb momentum force. During the first half second the 
calculations yield a large vaporization force (comparable to surface 
tension forces) and represent an intolerable sit\xation as far as the 
free-float time is concerned. However, due to the large temperature 
gradients existing, the resulting surface tension driven convection 
will tend to reduce these values. 

The Skylab orbit effects are included in the analysis by the 
method of Parker and Gatev/ood [.9] . The computer program pre- 
sented in Reference 9 was used as a subroutine in the present 
analysis to couple the effects of M512 position and orbit to the physi- 
cal phenomena occurring in the vacuum chamber. 


E. CONCLUSIONS BASED ON TRAJECTORY COMPUTATIONS 

It has been shown that free-float times from 32 to 48 seconds 
(depending on the material; are required for complete containerless 
solidification to occur. The magnitude of each of the physical forces 
was discussed and general comments made concerning the effect of 
each on free-float time in the previous section. This subsection 
svimmarizes some of these conclusions. 

• Electron beam force for a' 460-g spring load has 
a duration of approximately 100 milliseconds and 
gives rise to a float time of 40 seconds. 

• Kinetic energy of sphere due to spring effects 
(deployment technique) gives rise to float times 
of gr' ater than 50 seconds. 

• Vaporization force produces float times that a: 
probably the critical link in process. Best esti- 
mates of the corresponding float time between l/4 
second after eb cutoff and occurrence of an iso- 
thermal outer layer range from f to 20 seconds. 

For some cases the time could even be consider- 
ably less than 5 seconds. 

• Skylab orbit effects give rise to 29 to 46 seconds 
of float time. 
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Tablp 12 

SUMMARY OF VAPORIZATION FORCE CALCULATION 


Nil kol 


Vaporir.alion Tiino (sec) 

Force 1 

I 

1 


(flyin' s) 


3.00 





1 


1 


1 


Superposition of all of these forces produces an 
additive (as opposed to canceling) effect in most 
cases. Trajectory calculations show that even 
with conservative vaporization forces float times 
of less than 15 seconds will occur. 

I Lockheed documented in December 1972 10 an early projec- 

tion of anticipated float times less than the required containerless 
j solidification times and again in April via the Science & Engineering 

I Information Technical Note 11 , presented as Appendix A for reference, 

j In the latter report a conservative e. timate of less than 20 seconds was 

I forecast. 

' Anothei potential problem area, discussed in the 27 March 1973 

presentation at MSEC, is that even with the relatively large forces 
J on the sphere during the critical time oi deployment, the possibility 

of "sticking" exists. Films of the KC-135 flight test showed that an 
adherence to the ceramic occurs on occasion that finally results in 
' deployment due to acceleration transients in the ballistics trajectory 

of the aircraft. A typical sequ r :e is shown in Figure 46. The events 
are referenced to i.ne time corresponding to what appeared to be the fully 
molten condition, t=0. For the ensuing 0.300-secc ■’d, that adherence 
to the ceramic occurred with periodic motion due to the acceleration 
transients is evident in the film sequence. At 0.310-second, deploy- 
ment occurred with the trajectory shown in Figure 46. The important 
point remains, however, that without the acceleration transients 
(relatively large compared to the Sky lab gravitational environment) 
specimen retention could have occurred. This adhesion is probably 
due to contamination of the ceramic from previous melts. The con- 
clusion is that the possibility of unintentional retained specimen may 
occur in the actual Skylab mission. 



F . LIQUID DYNAMICS 

This subsection presents a summary of calculations made for 
Ni and Ni-12% Sn based on the thoughts outlined in Section II. First, 
the formation times are shown in Table 13 where it is seen that the 
transformation to minimum energy occurs very rapidly. However, 
due to inertial considerations, the spherical free surface is "overshot" 
and oscillations are set up which continue until dissipated by viscous 
damping. Table 15 gives the oscillation frequencies while Table 14 
shows the time required for viscous effects to dampen an initial 
perturbation to 0.01% and 1% of its initial value. 
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Table 13 

FORMATION TIMES FOR M553 MATERIALS 


M553 Materials 

P 

a 

R 

T 


(g/ cm^) 

(dyne/cm) 

(cm) 

(seel 

Ni 

7.85 

2050 

.3175 

.011 

Ni-12% Sii 

7.80 

1500 

.3175 

.013 


Table 14 

FREQUENCY FOR M553 MATERIALS 


Material 

Frequency 

(Hz) 

Ni 

41 

Ni-12% SN 

35 


Table 15 

DECAY TIME FOR M553 MATERIALS 


Material 

P 

P 

R 

Time 

Time 

IB 


(g/cm-sec) 

(g/cm^) 

(cm) 

(constant) 

(.01%) 


Ni 

.05 

7.85 

.3175 

3.14 

28.8 

14.5 

Ni-12% Sn 

.043 

7.80 

.3175 

3.62 

33.0 

16.7 






























t = .295 (sec) t = .310 (sec) 

FIGURE 46. TRAJECTORY .SEQUENCE - MARCH 1972 KC-135 FLIGHT 
(NICKEL) (CONTINUED) 


93 







FIGURE 46. TRAJECTORY SEQUENCE (CONCLUDED) 
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G. SOLIDIFICATION PATTERNS 


The theoretical preferential solidification pattern can be de- 
termined by determining the order that the individual nodes reach 
the solidi^ic^tion temperature on the cooldown cycle. For deployed 
specimens ii has been shown that isothermal layers develop within 
one second after eb cutoff and as shown in Figure 31. Solidification 
occurs from this outer layer inward as expected. The time that 
solidification begins on the surface is seen to be 10.5 seconds com- 
pared to 10.7 seconds for the center of the sphere as there are only 
a few degrees difference between the outer and inner portion from 
4 to 10 seconds after eb cutoff. However, the retained specimen 
exhibits an entirely different phenomenon. The initial solidification 
point is the sting as expected. The pattern develops (namely, solidi- 
fying from the outside in) as areas 1-2-8, 3-4-8 and 5-6-8 begin to 
solidify. As areas 7 and 8 reach the solidification temperature, 
area 9 is still superheated; area 11 is, however, molten at this point 
and no shrinkage should occur. The same is true of area 10-11-8, 
but as area 12-13-8 approaches the solidification temperature at 
area 13, shrinkage could occur as this is the last region to begin the 
solidification process — thus a possible void. This occurs as a result 
of the relatively large vaporization losses in the vicinity of the eb 
impingement, giving rise to large cooling rates. 

These calculations were made with the eb impingement location 
corresponding to actual operating procedures. This corresponds to 
only one-half of the beam impinging on the surface in the region of 
area 10-12 in Figure 47. Additional calculations were made with full 
eb impingement in the center of the specimen 7-5. It was found that, 
depending on the power (or beam duration after sting release), the 
region opposite the heat addition point (lower left of specimen in 
Figure ^7) would not melt before deployment. Then depending on the 
amount of superheat the remaining solid could melt in the free-float 
condition. Alternately, if eb impingement was lower than the center 
and sting release occurred early in the melt cycle, a large unmelted 
portion would remain in the final product. As stated previously if 
eb impingement is near the "one -o' clock” position, a fully molten 
specimen will result. 


i 

i • 1 






Deployed 


Retained 


FIGURE 47. SPHERE SOLIDIFICATION SEQUENCE (PURE NICKEL) 






SECTION V. CONCLUSIONS AND RECOMMENDATIONS 


'i Analyses of both flight and ground-based tests of the M551 

metals melting experiment indicated no definite trends in the effects 
of gravity on the test results. The beading phenomenon which appears 
^ in the ‘Stainless steel and tantalum specimens is shown to exhibit a 

remarkable similarity to vortex shedding. The bead spacing distance 
predicted from vortex shedding theory was found to be in almost exact 
agreement with test results for the stainless steel specimens and within 
about a half order of magnitude for the tantalum specimens. 

In the M553 sphere forming experiment, the predicted tempera- 
ture histories compared very closely to ground-based test results. 

This experiment proved to be more difficult to achieve satisfactory 
results than the M551 experiment. Some of the sphere specimens 
did not deploy from the sting, and the intense localized heating and 
pressure from the electron beam induced non-uniform temperature 
distributions and surface oscillations. Also, the deployment mech- 
anism and the pressure forces induced by the electron beam resulted 
in sizeable drift velocities. It is suggested that future improvements 
in experiment design for this type of study might include indiiction 
heating to provide a more uniform temperatiire distribution and to 
eliminate surface deforming pressures, and magnetic confinement to 
eliminate drifting of the specimen into the chamber wall. 
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FOREWORD 

This technical note was prepared by the Fluid Mechanic Applications Group, Aerd- 
I Mechanics Section of the Lockheed-Huntsville Research & Engineering Canter for \ 

the Process Engineering Laboratory of Marshall Space Flight Center under Contract 
NAS8-Z87Z9. This task was prepared at the request of Mr. Paul H. Schuerer, ! 

S&E-PE-MX. 

INTRODUCTION 


The information contained herein is for the purpose of updating the documentation 
for Skylab experiments M551 and M553. To incorporate the latest observations 
and opinions of contractors supporting development of the experiments and evalua- 
tion of resulting specimens and data, the following questions are asked: 



1. What do you consider to be the most significant results of your 
studies to date? 

Z. What are your projections for variations -in-sample properties, 
composition, structure and experiment data attributable to 
"zero gravity effects"? 

3. What specific properties or effects will you be concentrating 
your efforts on in evaluating the flight samples and data? 

4. How can the proposed variations noted in (Z) above be used 
to structure additional experiments, improve commercial 
materials and processes and develop new materials? List 
specific examples. 

5. What particular statements should be included in the experiment 
objectives as a result of your specific studies and projections 
of zero-g effects? 


Each question is answered for both the M551 metals melting and the M553 sphere 
forming experiments. 
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DISCUSSION 
Question 1 

Gene rai : Lockheed -Huntsville is responsible for all process phenomena with 
particular attention given to adhesion -cohesion studies with emnh'asis v-> the 
fluid dynamics of the molten metal. The effects of cohesion of vhe molten metals 
and their capability of relating their shapes have been interrelated with the 
adhesive influence of the molten metals in contact with ceramics and metals in 
a low -gravity environment. 

Also the importance of fluid dynamic considerations has been established in 
that all expected variations in terrestrial versus space processing are created 
as a result of the dominance of capillary motion (instead of gravity-induced 
flow), minimal segregation (sedimentation) and a decrease in hydrostatic pres- 
sure gradient — all hydrodynamic phenomena. This conclusion can be restated 
that gravity has no direct effect on grain structure or any other property of the 
solidified material. The gravity effects all occur as a result of differing fluid 
motion. 

• M551 

e Establishment of the role of vaporization in the cavity shape 

e Performance of three-dimensional temperature calculations for 
weld and dwell modes of melting for multi -component variable 
property material 

e Establishment of magnitude and direction of gravity during actual 
Skylab mission to allow meaningful conclusions to be drawn from 
micrographs about zero-gravity effects. 


• M553 

e Identification of physical forces affecting the trajectory of the 
molten sphere 

e Development of unifying calculational technique to predict tra- 
jectory of molten spheres during free -float condition 

e Development of three-dimensional thermal model to predict 
temperature history of sphere during melting 
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• Establishment of the role of vaporizatior on temperature history 
and motion of sphere 

e Establishment of the magnitude and qualitative motion of surface 
dynamics 

e Performance of detailed analysis of KC-135 flight data tmd ground- j 

based experimental tests. 

t 

Question 2 

General ; As discussed previously, the only expected differences in zero-g 
versus one-g processing are the dominance of capillary motion instead of 
gravity-induced flow, the absence of segregation and a decrease in hydrostatic 
pressure. There are, however, many indirect effects (advantages) due to a 
lack of gravity. In general, the factors affecting the final sample properties, 
c imposition and structure can be divided into three major categories. i.e.. 

e Nucleation 
e Growth 

e Reaction kinetics. 

The absense of gravity might allow better control of the indep - ^ ^ variables 

affecting each category: this occurs directly through the three atorementioned 
fluid dynamic effects and indirectly through other mechanisms that are a function 
of fluid motion (gas content, composititm, foreign particles, imperfections, cool- 
ing rate, nucleation rate, etc.) 

• M553 

• Less heterogeneous nucleation will occur due to foreign particles 
as segregation effects will be minimal. 

e If the material has any gas content that comes out of solution during 
solidification, an increase in porosity may result as there will be 
minimal tendency for the gas pockets to surface (except for convec- 
tion effects). 
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• For the materials with high viporUation rates (Ni> 12% Sn, imd 
Ni-30% Cu), a substantial increase in surface cooling occurs 
thus greater nucleation and surface tension gradients. 

• Even though the cooling rate is dependent on the amount of con- 
vection, the cooling rate will be determined by radiation con- 
siderations only as the sphere reaches isothermal superheat 
conditions in approximately a second. This in turn dictates the 
solidification rate* 

• Since solid formatior in poly crystalline materials is highly 
dependent on the distribution of foreign particles and the liquid 
motion, the associated differences should be discernable (homo- 
geneous nucleation). 


• M551 

e Since all fluid mechanic phenomena require a finite time to occur 
(due to inertia considerations) an effect due to gravity during the 
weld portion of the experiment may not occur (except for a more 
spherical bead when beading occurs) as large cooling rates are 
available via conduction. 

e A more homogeneous dwell structure will occur due to lack of 
sedimentation. 

e The liquid cavity shape will be determined by vaporization effects 
and should be nearly symmetric (different temperature gradients 
do exist in positive and negative radial directions). 

e If degassing is the cause for the splattering as noted in color 
films of the process, porosity may result in the zero-gravity 
dwell as the driving force for the pockets to surface is less. 

Question 3 

• M551 and M553 


e The flight films will be analyzed in detail and compared with KC- 
135 and ground-based results. 

e Observations from other contractors will be analyzed and compared 
with expected and predicted phenomena. 

• Correlation of actual gravity direcvion during Skylab experiments 
with experimental observations and results. 
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Question 4 
• M553 

To achieve the objectives of the current M553 experiment for all the selected 
materials, free -float times on the order of 60 seconds are required. Since it 
is highly unlikely that all of the physical forces will add in such a fashion to 
negate each other (as the system was not designed to do so), actual float times 
will be much less (probably less than 20 seconds). However, the concept of the 
sphere forming experiment is extremely exciting as the production of a nearly 
perfect homogeneously structured spherical specimen has important terrestrial 
applications. 

• Future experiment design should emphasize non -mechanical 
deployment, uniform heating and cooling aind filming rates 
of at least 3000 pps. 


O M551 

• The results of the weld portion of the experiment will establish 
the feasibility of "space repair." 

• Additional experiments with the same unit with high-speed color 
mo\n( s (3000 pps) with a different lens system should allow the 
fluid mechanics of welding to be analyzed in sufficient detail to 
answer mainy currently unresolved problems. 

Question 5 


• M553 


• Emphasis should be placed on fluid mechanic consideration as 
outlined previously. 

• Study of deployment dynamics as experiment success and final 
product are both dependent on this important phase of the ex- 
pe riment. 

• Role of temperature history on the final product. 
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• M551 


• Role of temperature history on the final product. 


" '/■ j ( Hi J 

M. R. Drashears, Ph.D. 

Fluid Mechanic Applications Group 


Approved 

B. Hobson Shirley, Supervisor 
Aero 'Mechanics Section 


A-6 


I 



1 


1 


1 


I 


Appendix B 

THERMAL HISTORY RESULTS FOR ALUMINUM 
AND TANTALUM DISKS 


> 


I 


i 



Appendix B 


Temperature calculations were made for the aluminum and tantalum 
disks in the same manner as for the stainless steek disks described in 
Section 3. The disk configuration (thickness, etc.), for all disks was as- 
sumed to be identical. The resulting temperature plots are given in Figs. 
B-1 through B-20. The maximum calculated temperatures for all three 
disk materials are compared in Table B-1. 


Table B-1 


COMPARISON OF CALCULATED MAXIMUM TEMPERATURES 
FOR STAINLESS STEEL, ALUMINUM AND TANTALUM DISKS 


Node 

No. 

Maximum Temperatures 

°C 

Stainless 

Steel 

Aluminum 

Tantalum 

317 

322 

697 

879 

315 

236 

669 

880 

223 

330 

198 

126 

221 

332 

201 

126 

219 

37 

220 

117 

135 

725 

1010 

1061 

133 

725 

1092 

1060 

131 - 

41 

218 

152 

35 

818 

1256 


33 

819 

1278 


31 

58 

219 

157 
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FIGURE B-l. THERMAL HISTORY FOR A 0,02-tNCH ALUMINUM 
DISK (NODES 4, 10, 16, 28, 40 AND 52) 
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FIGURE B-2. THERMAL HISTORY FOR 0.02-INCH ALUMINUM 
DISK (NODES 31 THROUGH 36) 
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FIGURE B-3. THERMAL HISTORY FOR 0.05-INCH ALUMINUM 
DISK (NODES 104, 110, I 16, 128, 140 AND 152) 
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FfGURE B-4. THERMAL HISTORY FOR 0.05-INCH ALUMINUM* 
DISK {NODES 131 THROUGH 136) 
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FIGURE B-5. THERMAL HISTORY FOR 0.125-INCH ALUMINUM 
DISK (NODES 219 THROUGH 224) 
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FIGURE B-6. THERMAL HISTORY FOR 0.125-INCH ALUMINUM 
DISK (NODES 249 THROUGH 254) 
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FIGURE B-7. THERMAL HISTORY FOR 0.125-INCH ALUMINUM 
DISK (NODES 279 THROUGH 284) 
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FIGURE B-8. THERMAL HISTORY FOR 0.125-INCH ALUMINUM 
DISK (NODES 301 THROUGH 306) 
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FIGURE. B-9. THERMAL HISTORY FOR 0.250-INCH ALUMINUM 
DISK (NODES 315, 345, 375, 317. 347 AND 377) 
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FIGURE B-10. THERMAL HISTORY FOR 0.250-INCH ALUMINUM 
DISK (NODES 310, 322, 328, 316, 346 AND 376) 
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FIGURE B-ll. THERMAL HISTORY FOR 0.02-INCH TANTALUM 
DISK (NODES 4, 10. 16, 28, 40 AND 52) 
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FIGURED- 12. THERMAL HISTORY FOR 0.02-INCH TANTALUM 
DISK (NODES 31 THROUGH 36) 
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FIGURE B-13. THERMAL HISTORY FOR 0.05-INCH TANTALUM 
DISK (NODES 104, 110, 116, 128, 140 AND 152) 
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FIGURE B-14, THERMAE HISTORY FOR 0.05-INCH TANTALUM 
DISK (NODES 131 THROUGH 136 ) 
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FIGURE B-15. THERMAL HISTORY FOR 0.125-INCH TANTALUM 
DISK (NODES 219 THROUGH 224) 
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FIGURE B-16. THERMAL HISTORY FOR 0.125-INCH TANTALUM 
DISK (NODES 249 THROUGH 254) 
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FIGURE B-17. THERMAL HISTORY FOR 0.125-INCH TANTALUM 
DISK (NODES 279 THROUGH 284) 
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FIGURE B-18. THERMAL HISTORY FOR 0.250-INCH TANTALUM 
DISK (NODES 301 THROUGH 306) 
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FIGURE B-19. THERMAL HISTORY FOR 0.250-INCH TANTALUM 
DISK (NODES 315, 345, 375, 317, 347 AND 377) 










